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ABSTRACT 

We are looking at quasar environment in the context of large-scale structure (LSS) — 
a new approach, giving a more informed interpretation of quasar-galaxy correlations. 
This paper presents our first results for a sample of 2; ~ 0.3 quasars. We are using 
Voronoi tessellation applied in colour {Bj — R) slices for the detection of galaxy clusters 
and the minimal spanning tree (MST) to delineate the large-scale structure. This new 
cluster-detection method allows us to find reliably galaxy clusters at z < 0.3 from 
SuperCOSMOS measurements of UK Schmidt plates. By reconstructing the large- 
scale structure in a relatively narrow redshift band (0.2 < z < 0.3), we show that 
quasars follow the large-scale structure traced by galaxy clusters. None of the quasars 
in our radio-quiet sample is located in the central area of a galaxy cluster. Two quasars, 
found in a very rich environment, are actually located between two very close galaxy 
clusters, consistent with results on z ~ 1 quasars suggesting that cluster mergers may 
be involved in one of the quasar formation mechanisms. 

Key words: galaxies: clusters: general - quasars: general - cosmology: large-scale 
structure of Universe. 



1 INTRODUCTION 



Quasar environments have been studied on different scales, 
from those of the host galaxy to those of the quasar-cluster 
cross-correlation function. Yee & Green (1984) found that 
there is a difference in richness between the environments of 
radio-quiet and radio-loud quasars (RQQs and RLQs respec- 
tively). This result has been confirmed by EUingson, Yee & 
Green (1991), showing that RQQs are much less frequently 
situated in very rich galaxy clusters than RLQs. Confirma- 
tion came also from Hintzen, Romanishin & Valdes (1991), 
who found a significant excess of galaxies in the vicinity of 
RLQs, and from Boyle & Couch (1993), who found no sig- 
nificant excess of faint galaxies in association with RQQs. 
Hutchings, Crampton & Johnson (1995) obtained a differ- 
ent result; they observed the galaxy environment of RLQs 
and RQQs and concluded that there is no significant differ- 
ence in the richness. Since the RQQs in that sample were 
members of a large quasar group (LQG) (Crampton, Cow- 
ley, & Hartwick 1989, 1987), this result raised the question of 
whether the environments of members of LQGs could be dif- 
ferent. However, some of the more recent results for quasars 
in general imply that the richness of the environments of 
radio- loud and radio-quiet quasars could indeed be very sim- 
ilar (e.g. Saxton, Hall & Turner 1999, Wold et al. 2001), 
whereas others still support the scenario of RLQs inhabit- 



ing more dense groups of galaxies than RQQs (Hutchings et 
al. 1999, Teplitz et al. 1999, Sanchez & Gomzales-Serrano 
1999). 

The above results have been obtained using galaxy 
counts or the quasar-gala:xy cross-correlation function 
within small radii (< lh~^ Mpc) around the quasars. The 
interpretation of the results is not a straightforward task, 
because of the erasure of directional information. For exam- 
ple, quasars residing centrally in poor clusters would pro- 
duce similar statistics to those on the priepheries of very 
rich clusters. 

Originally, the twofold environment of quasars was in- 
terpreted as supporting the connection of RLQs with ellipti- 
cal host galaxies and RQQs with spirals (e.g. EUingson et al. 
1991). However, in the past few years evidence has emerged 
that RQQs may have diverse host galaxies. In a sample stud- 
ied by Bahcall et al. (1997), for example, more RQQs have 
been found in elliptical hosts than in spirals. The result has 
been confirmed by Boyce et al. (1998), and by McLure et al. 
(1999), who found that, just like radio-loud quasars, essen- 
tially all RQQs with luminosities greater than Mji = —24 
reside in massive ellipticals. In many cases the host galax- 
ies show traces of continuing interactions or past mergers 
(e.g. Bahcall et al. 1997, Boyce et al. 1998, Hutchings et al. 
1999). Even if the host galaxies appear in the optical pass- 
bands to be isolated and apparently undisturbed, they are 
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found to exhibit continuing or remnant tidal HI disruptions, 
indicating galactic encounters or mergers (Lim & Ho 1999). 

The emerging picture is that RLQs have elliptical host 
galaxies and prefer a cluster environment whereas the RQQs 
inhabit various host galaxies in a wide range of environ- 
ments. At present we have no information on how the close 
environment (< lh~^ Mpc) of quasars is related to the large- 
scale structure traced by galaxy clusters, and if indeed the 
local enhancements of the density in the vicinity of quasars 
can be interpreted as galaxy clusters. This prompted us to 
introduce a new approach and look at the environments of 
quasars in the context of large-scale structure (LSS). 

At low redshift, the LSS can be conveniently recognised 
as density peaks in multidimensional parameter space (ra, 
dec, magnitude, colour) (Gladders & Yee 2000). Wide cov- 
erage (over 6x6 deg^) is required to ensure the isotropic 
detection of structures on scales ~ 60h~^ Mpc. In this pa- 
per we describe our first results from this new approach. 
The quasar and galaxy samples used in this study are de- 
scribed in Section 2. The galaxy cluster detection method, 
including a short introduction to Voronoi tessellation, is out- 
lined in Section 3. In Section 4 we have compiled the clus- 
ter catalogue from which the LSS has been reconstructed. 
The results of the study of the environments of quasars are 
contained in Sections 5 and 6. The advantages of the new 
approach and the implications of our results are discussed in 
Section 7. Throughout this paper we adopt a standard cos- 
mological model with Q = 1, and Ho = lOO/i km s~^ Mpc~^. 



2 DATA 

Our quasar sample consists of 15 low-redshift {z < 0.4) 
quasars from two surveys covering most of ESO/SERC field 
927 (1950 field centre lO'' 40"^ 00= +05° 00' 00"). Five 
quasars are from the AQD survey which uses Automated 
Quasar Detection (AQD — Clowes 1986; see also Clowes, 
Cooke & Beard 1984). The spectroscopic observation of 
high-grade candidates were carried out at the CTIO 4-m 
telescope (Clowes & Campusano 1991; Clowes & Campu- 
sano 1994; Clowes, Campusano & Graham 1999). Eight 
quEisars come from the Chile-UK Quasar Survey (CUQS — 
Newman et al. 1998; Newman 1999) which uses the ultravi- 
olet excess (UVX) method for selection. The spectroscopic 
observations were carried out at the Las Campanas 2.5m Du 
Pont telescope. Two quasars have come from Keable (1987). 
The CTIO quasars have faint limits {Bj ~ 20.5 mag) and 
sparse sampling whereas the LCO quasars have bright limits 
{Bj ~ 19.5 mag) and complete sampling. The quasar sample 
will be incomplete, but, nevertheless, free of any spatial bias. 
Table 1 contains the basic information about the quasars in 
our sample. 

Objects for the determination of the density distribu- 
tion have been drawn from UKST Bj and R plates covering 
^ 6.3 X 6.3 deg^ of sky. The plates have been digitized on 
the SuperCOSMOS measuring machine at the Royal Obser- 
vatory Edinburgh. Only objects recognized on both plates 
became a part of the final catalogues. This process elimi- 
nates virtually all of the spurious detections in the individ- 
ual plate scans (due to satellite and meteor trails, emulsion 
noise, etc.). The blue plates (J10128 and J10063) have been 
calibrated by applying 3'''^ degree polynomial fitting to deep 



(limiting Bj ~ 21.5) CCD sequences from Graham (1997). 
For the faint end {Bj > 21.5) a linear extrapolation has 
been apphed. The R plates (R18775 and R10071) have been 
calibrated in a similar way but using Keable (1987) data as 
secondary calibrators. The photometrically calibrated data 
set is limited at Bj = 22.5 and R = 21.5. 

Star-galaxy separation is not reliable at such faint lim- 
its, but for the purpose of this study faint stars can be 
treated as forming a uniform background. Jones et al. (1991) 
found that the distribution of star counts in cells (small ar- 
eas of size 4.96 x 4.96 arcmin^), at similar galactic latitude 
to our study {I ~ 50°), is close to that expected for a homo- 
geneous Poisson distribution and that the variations of star 
counts are far lower than those of the galaxy counts. Consid- 
ering that there are about three times as many galaxies than 
stars at the limiting magnitude {Bj = 22.5) the assumption 
about the uniformity of the star distribution should be a 
good approximation. 



3 FINDING GALAXY CLUSTERS 

Galaxy clusters arc the most massive gravitationally bound 
objects known, and are the obvious tracers of the LSS. They 
appear as peaks in the number density distribution of galax- 
ies. Since we are interested in the position of quasars with 
respect to cores of galaxy clusters and want to resolve the 
subclusters if present, it is crucial to preserve the topo- 
logical information when sampling the density distribution. 
The Voronoi tessellation has been chosen, since it is a non- 
parametric method. Voronoi tessellation, also known in com- 
putational geometry as Dirichlet tessellation, provides a par- 
tition of a point pattern according to its spatial structure. 
Features of this kind of decomposition can also be used for 
analysis of the underlying point process and provide a non- 
parametric approach to complex point processes, such as 
galaxy clustering. 

Given a set S of n distinct points in TZ"^, the Voronoi 
diagram is the partition of??.'* into n polyhedral regions a{p) 
for (p £ S). Each region of a{p), called the Voronoi cell of 
p, is defined as the set of points in TZ"^ which are closer to p 
than any other points in S: 

a{p) = {xen'^ I f{x,p) < f{x, q) V q € S - p} 

where f is the Euclidean distance function (Okabc et al. 
2000 and references there). One can use different distance 
functions to define various variations of Voronoi diagrams, 
but they will not be discussed here. The set of all Voronoi 
cells and their faces forms a cell complex. The vertices of 
this complex are called the Voronoi vertices, and the ex- 
treme rays (i.e. unbounded edges) are the Voronoi edges. 

A useful application of the Voronoi model is the compi- 
lation of surface density maps from point data. If one gen- 
crates the Voronoi diagram of the point set S and measures 
the area Ai of each polygon a^, one can consider as 
an indicator of the local intensity of the point pattern at 
ai (Duyckaerts, Godefroy & Hauw 1994). Using the astro- 
physical nomenclature, Ai is the local number density at the 
position of the object i from the sample. 

The Voronoi tessellation diagrams permit the generar 
tion of contour plots of the density distribution without 
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Table 1. Low redshift quasar sample. X and Y are the coordinates of the quasars in arcmin- 
utes from the the field centre. The notation AQD/UVX means AQD quasars re-discovered 

by CUQS. 



RA (J2000) DEC (J2000) X Y z Bj Survey 

arcmin arcmin mag 



10 


33 


56.8 


+05 


48 


30 


-130 


64 


0.376 


19.3 


AQD/UVX 


10 


34 


47.9 


+05 


16 


16 


-117 


32 


0.331 


19.2 


AQD/UVX 


10 


34 


58.3 


+05 


52 


32 


— 114 


68 


0.297 


19.5 


uvx 


10 


35 


06.0 


+06 


01 


41 


-112 


77 


0.245 


18.9 


uvx 


10 


36 


15.7 


+03 


19 


15 


-95 


-85 


0.389 


18.6 


AQD 


10 


36 


26.3 


+04 


54 


35 


-92 


10 


0.382 


18.3 


AQD 


10 


38 


28.2 


+04 


45 


35 


-62 


1 


0.286 


18.9 


uvx 


10 


38 


31.1 


+06 


09 


12 


-61 


85 


0.290 


18.1 


uvx 


10 


38 


57.9 


+04 


48 


39 


-54 


1 


0.357 


19.0 


uvx 


10 


44 


06.4 


+03 


14 


50 


23 


-89 


0.267 




Kcable 


10 


44 


52.0 


+03 


52 


51 


34 


-51 


0.213 


17.2 


Keable 


10 


45 


04.8 


+04 


34 


01 


37 


-10 


0.395 


19.5 


uvx 


10 


45 


35.1 


+03 


41 


22 


45 


-63 


0.355 


19.3 


UVX 


10 


46 


10.6 


+03 


50 


31 


54 


-54 


0.363 


19.2 


UVX 


10 


50 


42.8 


+04 


17 


39 


122 


-27 


0.230 


19.4 


AQD/UVX 



gridding the data. Gridding of data introduces, in many 
cases, undesired smoothing and parameterisation, which dis- 
tort the morphology of the measured distribution. 

The Voronoi tessellation has been used successfully for 
finding galaxy clusters (Ramella et al. 2001, Kim et al. 2000) 
and its performance compared with that of the matched fil- 
ter method by Ramella et al. (2001). In contrast to Ramella 
et al. (2001), we detect galaxy clusters by thresholding of the 
density contours in colour slices and not in magnitude slices. 
The use of colour slices has been proposed and tested by 
Gladders & Yee (2000). Selecting galaxies in colour slices in 
the colour-magrutude diagram increases the density contrast 
of any grouping of eaxly-type galaxies, which are normally 
found in clusters or groups. A colour slice in the colour- 
magnitude plane provides separation in redshift space via 
the red sequence of early-type galaxies in clusters (Gladders 
& Yee 2000, and references inside), thus reducing consider- 
ably the projection effects. 

The redshift range < z < 0.4 corresponds to a Bj — R 
colour space in the range 1.4 < Bj — R < 2.5 using the 
galaxy colours derived by Fukugita, Shimasaku & Ichikawa 
(1995). The limit in the redshift space for this pair of filters is 
z ~ 0.4, since the colours for higher redshift ellipticals start 
to drop for z > 0.6, reaching Bj — R ^ 2.5 again at 2; ~ 
0.8. Such high redshift clusters are in any case beyond the 
magnitude limit of our data. Using the colours of the known 
clusters and the model predictions (Kodama et al. 1998) the 
appropriate gradient of the slices in the colour-magnitude 
plane has been deduced to be —0.06. The width of the slices 
is mainly dictated by the uncertainities in the photometric 
calibration and empirically selected to be A{Bj — R) = 0.2 
mag. The colour slices are overlapping in steps of 0.1 mag. 

The analysis is performed initially in every colour slice 
separately. The points belonging to candidates for galaxy 
clusters are selected as those exceeding an intensity thresh- 
old -I X^r^i where Ui is the area of the i*^ polygon 
(corresponding to the i"* object) and n is the total num- 
ber of points (polygons). This threshold corresponds to 
{N field + Nciuster)/N field ~ 2 where Nfieid is the number 
density of field galaxies integrated in the line of sight to 



2000 



1 500 - 



I 1000 - 

T3 




density [a rcmi n ^] 

Figure 1. The histogram illustrates the distribution of the 
Voronoi densities in a typical colour slice (in this example 1.9 < 
{Bj — R) < 2.1). The dashed vertical line indicates the threshold 
density. 



z = 0.3 and Nduster is the number density of the core of 

a poor cluster. The literature values for the densities have 
been adopted from Cox (1999). Figure 1 illustrates the dis- 
tribution of the densities in a typical colour slice, and the 
selected threshold, which usually picks ~ 5% of the cells in 
the slice. 

The colour-magnitude relation applies to the elliptical 
galaxies which in turn are known to form the cores of the 
clusters. The cluster cores can be expected to be mathemati- 
cally closed (no holes or single-object dips in the boundary) 
structures. Strongly non-virialised clusters, such as Virgo, 
consisting of multiple subclusters, may be considered as as- 
sociations of multiple separated cores. The algorithm em- 
ployed to close the boundaries of the cores is based on the 
principle that detected structures are closed if all cells with 
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cluster at X = 62.0423, Y = -99.9983 
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Figure 2. Two examples of detected clusters. For each cluster, the panel on the left is a contour plot of the density distribution, with 

the polygon enclosing the core of the cluster with early- type galaxies (triangles). The panel on the right is a colour-magnitude diagram 
of all objects contained in the polygon, with those objects indicated by asterisks defining the cluster detection. 



a simple majority of its natural neighbours (adjacent cells) 
exceeding the density threshold are being included in the 
selection. 

The Voronoi tessellation has the power to detect group- 
ings on the richness scale from loose groups to rich clusters. 
To find only groupings which qualify as clusters, a lower limit 
on the number of interconnected objects in a cluster must 
be introduced. The lower limit for a poor cluster (Abell class 
0) is by definition 30 galaxies. To detect such poor clusters 
with an ellipticals content as low as 25 per cent a minimum 
of 7 interconnected objects is required. The results from all 
the colour slices have been combined according to the prin- 
ciple that two clusters from different slices sharing two or 
more members constitute the same cluster. 



4 CLUSTER SAMPLE 

The application of the above procedure to the ~ 36 deg^ area 
of investigation produced a sample of 72 galaxy clusters with 
estimated redshifts in the range 0.03 < z < 0.37. The result- 
ing density is 2.0 clusters dcg^'^ with (z) = 0.19. Because 
of the magnitude limit, only 13 (optically brightest) clusters 
with z > 0.3 have been detected. Figure 2 shows some ex- 
amples of the detected clusters with different morphologies 
and redshifts. 

The data for all the detected clusters have been com- 
piled in Table 2. The positions axe expressed in RA and Dec 
(J2000) and as X and Y which are the coordinates in arc 
minutes from the field centre (used in all figures). The area 
of a cluster core is here defined as the area of the convex 
hull enclosing all the Voronoi cells of the early-type galaxies 
found to be members of the cluster. It gives a good measure 
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Table 2. Cluster sample. 
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Table 2. Cluster sample. Continued 
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of the core of the cluster, preserving its morphology. The 
number of galaxies (No) in a cluster core is here defined as 
the background corrected number of objects found within 
the convex hull defining the boundary of the core. The ra- 
dius (Rp) is the minimum radius necessary to enclose the 
whole convex hull. The colours of the red sequences have 
been established by inspection of the colour-magnitude plots 
for every cluster. The redshift z has been estimated using 
the red sequence approach and the colours of elliptical galax- 
ies compiled by Fukugita ct al. (1995). The uncertainty of 
the estimate is ~ 0.05 due to the uncertainties in the photo- 
metric calibration. The Abell richness (Nr) of the clusters 
has been computed using the estimated redshifts to calcu- 
late the angular extent of the Abell radius {1.5h~^ Mpc). 
The background counts, used for background corrections, 
have been taken cis those of a region without any detected 
clusters. The identifications of the previously known clusters 
have been indicated in the final column. 

The area of investigation contains 16 published galaxy 
clusters, 3 of which have known redshifts. We identified 14 of 
these clusters including all those with known redshifts. The 
2 clusters we could not identify may be spurious since they 
were detected only as density peaks in 2D or have few el- 
lipticals and so are below our threshold in the colour slices. 
Figure 3 gives the comparison of the published centres of 
the clusters and those determined using colour slices. The 
difference is strongest for the lowest redshift cluster (X=81, 
Y=22, z = 0.03); most of the other clusters agree very well. 
The three known redshifts (0.03, 0.07, 0.15) allowed an im- 
provement of our redshift estimate at z < 0.2 by direct cal- 
ibration of the colour-redshift sequence. 



5 THE LARGE SCALE ENVIRONMENT OF 

QUASARS 

The purpose of our investigation is to establish whether 
quasars trace the same structures as the galaxy clusters. 
The quasar sample covers the redshift range 0.2 < z < 0.4. 
However, the cluster sample appears to be very incomplete 
a,t z > 0.3 forcing us to limit the main investigation to one 
redshift slice, 0.2 < z < 0.3. The area of the quasar survey 
is smaller than that populated by detected clusters. Only 
objects in the common area (~ 25 deg^) of both quasar and 





□ 






□ 




□ □ 




□ 




□ 


□ 


□ + ° 




□ 




+ 


□ 


□ 

□ 


□ 


□□ 


n □ 


□ 




□ ° 


□ 



Bo □ 



ARA [arcmin] 

Figure 3. Cluster sample. The squares are centred on the galaxy 
clusters detected by Voronoi tessellation in (Bj ~ R) colour slices. 
The overplotted plus signs represent the known galaxy clusters 
drawn from the CEDAG Catalogue of Clusters of Galaxies. 



cluster samples have been selected for the environmental 
study. The final sample comprises 17 galaxy clusters and 7 
quasars (Figure 4). 

The problem of connecting discrete objects into a con- 
tinuous structure can be approached using another graph 
theoretical method, the minimal spanning tree (MST). This 
method has been successfully tested for finding LQGs in 
quasar samples (Graham 1997). The minimal spanning tree 
(MST) is a geometric construct, originating in graph theory, 
which was introduced by Kruskal (1956) and Prim (1957). 

In our case the nodes of the MST are the galaxy clusters. 
For the purpose of this study we did not attempt to separate 
the MST to recover the galaxy superclusters, but preserved 
the continuous LSS. When looking at Figure 4, it can be seen 
that quasars populate the region delineated by the MST as 
LSS. In order to test the statistical significance of this result 
we resampled the quEisar distribution using a Monte Carlo 
method. As a measure of how close the quasars follow the 
LSS, we have used the sum of the distances between quasars 
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Figure 4. The relation between the LSS in galaxy clusters (poly- 
gons) delineated by the MST (dotted line) and quasars (asterisks) 
in the rcdshift slice 0.2 < z < 0.3. The axes correspond to X and 
Y in units of arcminutes from the centre of ESO/SERC field 927. 
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Figure 5. The result of the Monte Carlo resampling of the quasar 
distribution in relation to the MST vectors. The histogram illus- 
trates the distribution of X/"— l*"*' ^^^'"'^ ^'^'^ distance of 
the i*** quasar to its closest MST edge, in 100,000 simulation 
runs. The vertical dashed line indicates the value produced by 
the quasars in the real sample. 



and their closest MST edges X)r=i ^t- The result of 100,000 
simulations has been plotted as a histogram in Figure 5. 
The test rejects the null hypothesis that the quasars are 
distributed independently of the LSS in clusters at a level 
of significance of 0.009. 
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Figure 6. The environment of the quasar at z = 0.245. The figure 
shows the density distribution (filled contours) and the cores of 
two clusters (polygons) identified at a comparable redshift. 

6 SMALL SCALE ENVIRONMENT OF 
QUASARS 

From Figure 4, it can be seen that the quasars appear to 
avoid the central density peaks of the clusters, and reside 
instead in two types of environments. Most of them (5 out 
of 7) reside on the distant (3 - 5h~^ Mpc) peripheries of 
the galaxy clusters tracing the same LSS. This is an impor- 
tant result, suggesting that the topology of the density field 
in the quasar environments might be a crucial parameter 
connected with the triggering of quasar formation and the 
controlling of their evolution. Some of the quasars (2 out of 
7), however, are found between two very close clusters. Such 
an arrangement suggests that quasar formation might also 
be associated with cluster mergers. 

Figure 6 illustrates the density distribution and the two 
close clusters around a quasar at z = 0.245. The estimated 
redshifts of both galaxy clusters arc comparable to that of 
the quasar, suggesting that all three objects are physically 
associated. 

The second merger candidates are more difficult to 
verify. Both clusters found in the vicinity of a quasar at 
z — 0.286 (Figure 7) arc obscured by foreground clusters. 
They have clear red sequences in the colour-magnitude di- 
agram but the enclosing polygons and redshift estimates 
are nevertheless less reliable than in the first case outlined 
above. The redshifts of both clusters are estimated to lie 
in the range 0.25 < z < 0.30 which is compatible with the 
redshift of the quasar. 

This result suggests that quasars found in rich environ- 
ments are actually not in the centres of galaxy clusters but 
mark the regions of cluster / subcluster mergers by resid- 
ing between two (or more) density peaks. With only two 
quasars, however, this needs to be confirmed with a larger 
sample (work in progress). 
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Figure 7. The environment of the quasar at z = 0.286. The figure 
shows the density distribution (filled contours) and the cores of 
two clusters (polygons) identified at a comparable redshift. The 
density enhancement superimposed on the position of the quasar 
is due to projection and does not form a cluster core. 

7 DISCUSSION AND CONCLUSIONS 

The majority of the studies of the richness of the environ- 
ments of quasars have relied on counting galaxies within 
a quite small radius ~ 1 Mpc without consideration of the 
morphology of the distribution. Such an approach provides a 
fast means of tabulating and comparing the richness within 
large data sets. However, to ensure a correct interpretation 
of the results, we must determine what physical structures in 
the galaxy distribution are associated with quasars. We have 
studied the association of quasars with optically detected 
clusters of galaxies. By reconstructing the LSS in a relatively 
narrow redshift band (0.2 < z < 0.3), we have been able to 
show that quasars follow the large-scale structure traced by 
galaxy clusters. The quasar-cluster LSS investigated in this 
paper is on the scale of superclusters (~ 60/i~^ Mpc), which 
suggests that quasars could be used to trace the distribu- 
tion of superclusters. Oort, Arp & de Ruiter (1981) brought 
forward for the first time the hypothesis that quasars are 
located in superclusters using the statistics of close quasar 
pairs as evidence. It is normally expected that clusters of 
galaxies trace mass, and that the LSS of galaxies also traces 
mass. Our results show that quasars trace the mass (galax- 
ies) on supercluster scales for low redshifts z ^ 0.3. It means 
that the formation/fuelling process of quasars is coupled, at 
low redshift, to the large density perturbations, which mates 
quasars useful indicators of meiss at z ^ 0.3. We should not, 
however, simply extend this conclusion to higher redshifts 
without observational support because the relative impor- 
tance of different mechanisms for quasar formation will al- 
most certainly vary with redshift. 

Secondly, none of the quasars in our radio quiet sam- 
ple is located in the central area of a galaxy cluster. Apart 
from the two quasars found between merging clusters, all the 
others reside on the distant (3 - 5h~^ Mpc) peripheries of 
galaxy clusters. The peripheral position of RQQs in relation 



to galaxy clusters was first noticed by Yee & Green (1984) 
when assessing by eye the distribution of companion galax- 
ies in the vicinity of quasars. Yee (1992) noticed that z ~ 0.6 
quasars do not reside close to the cD galaxies marking the 
centres of clusters, but these results do not seem to have 
been pursued further. Hcishimoto & Oemler (2000) found 
that galaxy interactions are more hkely in galaxy groups 
and very poor clusters. They found also about the same rar 
tio of interactions in rich clusters as in the field, three times 
less then in groups and poor clusters. The fact that most 
quasars in our sample (5 out of 7) avoid strongly the rich 
environments of the central regions of galaxy clusters is com- 
patible with models connecting the majority of low z quasars 
with galaxy groups and suggesting that major galaxy-galaxy 
mergers arc the main, but not the only, mechanism of quasar 
formation at low z. 



Lastly, we find that quasars found in very rich environ- 
ment (2 out of 7) may actually be located between two very 
close galaxy clusters. Due to the relatively small number of 
quasars in our sample it is not possible to draw conclusions 
about the percentage of such objects in the overall quasar 
population. Two clusters separated by a distance < 2h~^ 
Mpc are expected to have entered already the pre-merger 
stage (Schindler 2001) in which strong gas compression ac- 
companied by heating occurs. The high velocities of galaxies 
found in merging clusters do not favour major galaxy-galaxy 
mergers as a quasar formation mechanism, but rather point 
towards galaxy harassment (Lake, Katz & Moore 1998). 
However, cluster / subcluster mergers produce compression 
of hot gas between two subclusters shortly before they col- 
lide and shock waves after the collision (Schindler 1999). 
Considering the complex nature and the high energetics of 
the gas compressions and displacements already during the 
pre-merger stage of two galaxy clusters, we should consider 
also a model in which a quasar is formed in a galaxy-gas in- 
teraction. It has been argued that strong star formation can 
be triggered by the infiucnce of dense ICM at the confiucnce 
of a cluster / subcluster merger (Evrard 1991). Our result 
for z ~ 0.3 quasars, combined with the fact that Haines et 
al. (2001) and Tanaka et al. (2000) found quasars on the 
peripheries of clusters at z ~ 1.3 and 1.1, with some being 
accompanied by a band of star-forming galaxies and some 
evidence again for cluster mergers, suggests that such a pro- 
cess could indeed be one of the quasar formation mechar 
nisms. 



To determine the relative contributions of the cluster- 
merger like and distant periphery like environments to for- 
mation of quasars, we are planning to study a large quEisar 
sample, spanning a wide redshift range. The inclusion of 
radio-loud quasars will allow us to quantify the results as 
functions of redshift and radio-loudness. Multiwavelength 
observations will be required to establish the role of the ICM 
in triggering the formation of quasars. The X-ray morphol- 
ogy of the merging clusters will give us information about 
the distribution of the hot gas. Imaging of the host galax- 
ies and their close companions will be used to find traces of 
recent mergers or interactions. 
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